Infinite randomness fixed point of the superconductor-metal quantum phase transition 
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We examine the influence of quenched disorder on the superconductor-metal transition, as de- 
scribed by a theory of overdamped Cooper pairs which repel each other. The self-consistent pairing 
eigenmodes of a quasi-one dimensional wire are determined numerically. Our results support the 
recent proposal by Hoyos et al. (Phys. Rev. Lett. 99, 230601 (2007)) that the transition is charac- 
terized by the same strong disorder fixed point describing the onset of ferromagnetism in the random 
quantum Ising chain in a transverse field. 



Numerous recent experiments [l|, 0, S 0, [H, Q have 
measured the electrical transport properties of quasi-one 
dimensional nanowires. While thicker wires have vanish- 
ing resistance in the low temperature (T) limit, thinner 
wires do not display superconductivity even at the lowest 
T . The superconducting wires display clear signatures of 
thermal phase fluctuations of the Cooper pair order pa- 
rameter, >F, at low T. Quantum fluctuations of the phase 
and amplitude of \& increase with decreasing wire thick- 
ness, leading to a transition to a non-superconducting 
state. 

Recent work 0, H| has proposed that these experi- 
ments should be described by a quantum superconductor- 
metal transition (SMT) in the pair-breaking universality 
class. Arguments based upon microscopic BCS theory 
were used to propose a model of \& fluctuations damped 
by decay into single electron excitations of the metal 
Ha E3, EL 0, EH- In this paper, we will present the 
results of a numerical study of the influence of quenched 
disorder on this model. The role played by disorder near 
the quantum SMT is of considerable interest, as disorder 
correlations are of infinite range in the imaginary time di- 
rection and can lead to unusual critical phenomena 14 1 . 



In a renormalization group (RG) analysis of this over- 
damped Cooper pair model with disorder, Hoyos et al. 
[IH have recently argued that the SMT is described by 
a strong disorder fixed point. The latter exhibits acti- 
vated dynamic scaling where the logarithm of characteris- 
tic frequencies of \f r fluctuations grows as a power of their 
characteristic length scale. They argued further that the 
strong disorder fixed point is in the same universality 
class as the one describing the onset of ferromagnetism 
in the quantum random transverse field Ising model (RT- 
FIM) in one spatial dimension. Many exact results were 
obtained by Fisher [ll| for this fixed point, and carry 
through to a large extent to the case of the SMT. Note 
that this is a non-trivial result, since the RTFIM con- 
tains no dissipation, and possesses a discrete symmetry. 
All our numerical results below confirm the remarkable 
predictions of Hoyos et al. , providing strong evidence for 
the applicability of their strong randomness RG. 

Our analysis was carried out on a lattice discretiza- 
tion of the disordered overdamped Cooper pair model 
of Ref. at T = 0. The degrees of freedom are %(r), 



which are complex functions of imaginary time, t, on the 
sites, j, of a one-dimensional chain with action 
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where *&j(u>) is the Fourier transform of ^,-(r), and the 
couplings in S are all random functions of j. The quartic 
coefficients Uj are all positive to ensure stability and re- 
pulsion between Cooper pairs. The dissipation into the 
metallic bath is represented by jj , which is also required 
to be positive by causality. Finally, we can choose a gauge 
such that Dj > 0. A more careful analysis and suitable 
rescalings [17| allow us to reduce the randomness to the 
spatial dependence of Dj (uniformly distributed on (0, 1]) 
and ctj (taken to be Gaussian), while setting Uj — u and 
7j = 1. At zero temperature, the SMT can be tuned 
by reducing the mean of the ctj distribution, a, while 
keeping its variance constant at 0.25 in units of 7 2 . 

Equivalently, we can also work in a lattice model of 
fluctuating phases with ^fj(r) = e l6j ^ of unit magnitude 
[JEH; this should have the same properties as S, but our 
analysis proceeds more conveniently by also allowing for 
magnitude fluctuations. 

While S is a suitable model for describing the influence 
of disorder on the fluctuating Cooper pair states, we also 
have to consider the effect of randomness on the single 
electron states. We have estimated such effects in the 
framework of weak-coupling BCS theory: at criticality, 
we find that on a scale parametrically smaller than the 
single electron localization length, the gain in conden- 
sation energy can offset the cost in elastic energy when 
order parameter fluctuations take advantage of random- 
ness in the ctj. This justifies our focus on the influence 
of disorder in a purely bosonic overdamped Cooper pair 
theory. Details will be provided later (see also Ref. [ill) . 

The RG analysis [l5| was carried out in a model with 
an A-component order parameter and it was found that 
flows had only an irrelevant dependence on the value of 
N |18j. Thus the exact critical properties can be ob- 
tained by studying the model in the large N limit. This 
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is equivalent to approximating S by the Gaussian action 
S ° = E/ % \Pi\*i ^ +1 l 2 + ^ + l w DI*il 8 ] ' ( 2 ) 
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where the rj are determined self-consistently by solving 



-(I^W| 2 ) So 
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We set u = 1 to reach a strong coupling regime and 
use fixed but random boundary conditions, similar to 
those employed in Ref. 0- Solving the innocuous look- 
ing Eq. ([3]) for a large number of disorder realizations 
and large system sizes was the primary time-consuming 
numerical step in obtaining the results of this study. Sim- 
ilar numerical \aige-N methods have been used previ- 
ously for disordered s yste ms with conventional (power 
law) dynamic scaling [Til , liol ] but the presence of acti- 
vated scaling leads to sluggish dynamics and the neces- 
sity to properly include spurious disorder configurations 
that, although exponentially rare, can make large con- 
tributions to thermodynamic properties. The numerical 
solution is facilitated via the implementation of a method 
which we have dubbed the solve- join- patch (SJP) proce- 
dure. We begin by generating a realization of disorder for 
L sites with L large. Near the critical point, character- 
ized by the condition that the correlation length £ ~ L, 
the direct iterative solution of Eq. ([3]) is computationally 
quite costly. This is a result of the fact that the eigen- 
modes of So begin to delocalize and have a characteristic 
energy scale that is exponentially small in the distance 
from criticality, requiring that the solutions Tj must be 
computed with exponentially increasing precision. 

To cope with this difficulty, the system consisting of L 
sites is broken up into a group of smaller sub-systems 
with boundary conditions adjusted to reflect their lo- 
cation in the larger chain. The sub-systems are solved, 
then joined together in groups of two. The grouped sub- 
systems are now close to satisfying Eq. ((3|) and they can 
be quickly brought into accordance by patching, which in- 
volves re-solving a mini-system around the joint consist- 
ing of a small number of sites. The joined and patched 
sub-system can now be easily solved and the SJP pro- 
cedure is iterated until a full solution to Eq. ([3|) is ob- 
tained for the complete chain of L sites. We have consid- 
ered up to 3000 realizations of disorder for system sizes 
L = 16,32,64, and 128. 

Fisher's remarkable solution of the RTFIM [r| in- 
cludes asymptotically exact results for the exponents 
and correlation functions at the infinite randomness fixed 
point, and many directly translate to the RG calculations 
by Hoyos et al. 15] for the dissipative model considered 
here. In particular, one expects activated dynamic scal- 
ing with ln(l/f2) ~ where f2 is a characteristic energy 
scale and tp — 1/2 is a tunneling exponent. This reflects 
the fact that at an infinite randomness fixed point, the 
dynamical critical exponent z is formally infinite. The 



RG approach defines a real space decimation procedure 
that either creates or destroys clusters or bonds as the 
energy scale is reduced. The typical moment of a surviv- 
ing cluster scales like \x ~ ln^(l/f2) at criticality, where 
<p = (1 + y/5)/2 ~ 1.62 is the golden mean. Average 
correlations are described by a correlation length which 
diverges as £ ~ H" with v = 2 and S a measure of 
the distance from criticality. From Ref. [l5|, 6 is expected 



to be proportional to \n(ri/r c ) where r c is some critical 
value. Our numerical study reveals that close to critical- 
ity this quantity is linearly related to the detuning of the 
average a from its quantum critical value, a c (which has 
yet to be determined) and it further demonstrates that 
correlations among the due to their self-consistency 
does not affect the strong randomness RG flow. 

The remainder of this paper will present a numerical 
confirmation of the results of Ref. [l5j by providing argu- 
ments for the presence of dynamically activated scaling 
at the quantum SMT, characterized by exponents v, ip 
and <j> taking on their RTFIM values. The evidence comes 
from an analysis of equal time correlations, energy gap 
statistics and dynamic susceptibilities in the weakly dis- 
ordered quantum Griffiths phase [211 ]. 

We begin by studying the disorder averaged equal-time 
correlation function C(x) = (^>*(t)^o(t)} s , which can 
be computed from the quadratic effective action Sq once 
the full set of solutions {r^} has been obtained. In the 
disordered phase, where 6 = a — a c > the asymptotic 
form of C(x) for the RTFIM has been predicted to de- 
scribe both exponential as well as stretched exponential 
decay in addition to power law behavior 16] 



C(x) 



exp [-(x/Q - (27^ 2 /4) 1 / 3 (x/£) 1 / 3 ] 
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If we use this expression to define the correlation length 
£, we can perform fits for each value of L and various a 
to extract £(L, a) as is seen in Fig.Q]for L = 64. We find 
remarkable agreement (solid lines) with Eq. ([4]) over six 
orders of magnitude for all system sizes considered. 

As mentioned above, the length scale which describes 
average correlations is expected to diverge like £ ~ \&\~ v 
as the critical point is approached. We have employed 
this result to perform a log-log fit to the finite size scaled 
correlation length (data extrapolated to L — > 00) as a 
function of S, as is shown in the inset of Fig. [TJ The 
value of a c was found from the mean of the critical ctj 
distribution which minimized the least square error of 
power law fits involving S — a — a c . This leads to a value 
of a c = —0.93(3) for the critical point and v = 1.9(2) 
for the correlation length exponent with the number in 
brackets indicating the uncertainty in the last digit com- 
puted from the fitting procedure. The obtained value of v 
is in accord with the value of 2 predicted for the RTFIM. 
The correlation length could also have been defined via 
the exponential tail of C(x) at large separations which 
yields compatible values for both a c and v. 




FIG. 1: (Color online) The equal-time disorder averaged cor- 
relation functions for L — 64 and five values of the mean of the 
aj distribution, a. The solid lines are fits to the asymptotic 
form described in Eq. ((4]) via £ and an overall scale parameter. 
The inset shows the result of a fit to the power law form of 
the finite size scaled correlation length leading to an estimate 
for the location of the critical point a c = —0.93(3) and the 
correlation length exponent v = 1.9(2). 
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FIG. 2: (Color online) The finite size scaled value of the 
disorder averaged logarithm of the minimum excitation en- 
ergy plotted against the distance from the critical point 
S. We observe divergence consistent with the scaling form 
|hTn| ~ and using the value of a c and v found above we 

determine ip = 0.53(6) from a log-log linear fit (inset). 



For each realization of disorder and each value of a 
we define the gap fi(L) to be the smallest excitation en- 
ergy in the system, which in general corresponds to the 
most delocalized mode of So- Rare disorder configura- 
tions cause clusters to behave as if they were much more 
critical than the global value of 5 would suggest. These 
clusters dominate the critical modes and exhibit abnor- 
mally small gaps that make large contributions to dis- 
order averages of In 57, leading to the highly anisotropic 
scaling relationship between space and time that is the 
hallmark of strong disorder fixed points. An analysis of 
the probability distribution for the logarithm of the en- 
ergy gap in the RTFIM was carried out by Young and 
Rieger [22j where they found cogent evidence for z — oo. 
We have performed a similar analysis here, with the same 
result. In addition we find that lnfi, as the minimal ex- 
citation energy, is naturally Gumbel distributed. If ac- 
tivated dynamic scaling is indeed present, the disorder 
averaged value of the logarithm of the gap should scale 
like |lnfi| ~ £^ ~ 5~ u ^ where we have used the scaling 
form of the correlation length. Such divergent behavior 
for the finite size scaled value of | ln.fi | is demonstrated in 
Fig. [2j The possibility of conventional scaling was con- 
sidered but ultimately excluded through the examination 
of the maximum likelihood estimator for a wide range of 
power law fits. Using the previously determined values of 
a c and v, the tunneling exponent can be extracted from 



a log-log linear fit of the average logarithmic spectrum 
as shown in the inset of Fig. [2j producing ip — 0.53(6) 
which is consistent with the RTFIM prediction of 1/2. 

To confirm full agreement with the universality class 
of the RTFIM, we must finally determine the value of 
the exponent <f> which controls the average moment, 
/i ~ | bawl* 1 , of a cluster fluctuating with frequency u>. 
This can be accomplished by investigating the imaginary 
part of the disorder averaged dynamical order parameter 
susceptibilities after they have been analytically contin- 
ued to real frequencies. We are interested in the average 
(k = 0) and local susceptibilities defined by 



ImxM =Im -J2 (^HfoM), 



itu — >uj-\-ie 



Im xi^(w) = Im (m(iu)^j(iu)) 



iu)—*uj-\-ie 



(5) 



(6) 



where (• • • )s indicates an average with respect to the 
large- N action in Eq. @ as well as a site average. Note 
that u> is now a real frequency, and we point out that our 
facile access to such dynamical quantities is one of the 
perquisites of the numerical approach we have taken. All 
frequencies are measured with respect to an ultra-violet 
cutoff A w which is required for convergence when com- 
puting the set of solutions to Eq. (J3]). Physically, one can 
argue that at criticality, the average cluster moment will 
be given by the ratio of the average to local susceptibility 
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due to the extra sum over sites in Eq. ([5]) . We thus define 



#(w) - j / \ 



|lnw|* JF(<^|mu;|) , (7) 



and expect that the scaling function T approaches a con- 
stant when the dimensionless variable J^llnajl <C 1. 
In the quantum disordered phase with 8 V ^ \ \nui\ ^> 1, 
a scaling analysis predicts F(x) ~ x 1_c ^ [23j and hence 
R ~ ln W | [3, |23|. In order to determine the 

value of </>, it is useful to consider a rescaled value of 
the susceptibility ratio R{8) = R{uj)/{5 v ^\ In uj\) which 
should be frequency independent according to the pre- 
dicted scaling form for R(ui) as uj — ► 0. We plot the fi- 
nite size scaled susceptibility ratio in Fig. [3] for the three 
smallest values of <5, and find confirmation of its |lnw| 
dependence. The inset of Fig. [3] confirms the frequency 
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FIG. 3: (Color online) The real frequency dependence of the 
finite size scaled value of the disorder averaged susceptibility 
ratio defined in Eq. © for three values of the S = a — a c . 
We observe the predicted In a; behavior. After a suitable 
rescaling described in the text we find that R does not depend 
on frequency as u) — * (inset), and a log-log linear fit gives 
the value of the cluster exponent to be <f> = 1.6(2). 

independence of R and by determining the best linear fit 
of lni? to In 6 for w < ICE 3 with vip = 1.0(1), we find 
a cluster exponent <j> — 1.6(2) which is very close to the 
predicted RTFIM value of (1 + \/5)/2. 

The results of the above analysis, as highlighted in 
FigS-HHSl provide compelling evidence for the applicabil- 
ity of the real space RG analysis of Hoyos et ah, and 
further reproduces a number of results of [161 ] to un- 
expected accuracy. This confirms that the considered 
model for overdamped repulsive Cooperon fluctuations 



in the presence of quenched disorder near a SMT ex- 
hibits dynamically activated scaling and is controlled 
by an infinite randomness fixed point in the same uni- 
versality class as the RTFIM. The transition is charac- 
terized by the numerically computed critical exponents 
(i/,ip,<f>) ~ (1.9,0.53,1.6) which are entirely consistent 
with those of the one dimensional random quantum Ising 
model in a transverse field (2, 1/2, (1 + y/5)/2). 

In closing, we note that while our discussion has been 
framed in the context of the SMT, models similar to S 
describe the onset of a wide variety of orders in metal- 
lic systems [14(. Furthermore, the flow to the strong- 
disorder RTFIM fixed point is expected to also hold in 
higher dimensions [l5| . We thus propose that our results 
provide strong support for the applicability of the RT- 
FIM physics to many experiments involving the onset of 
spin- and charge-density wave orders in metals. 
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